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Maximizing the precision in estimating parameters in a quantum system subject to instrumentation constraints 
is cast as a convex optimization problem. We account for prior knowledge about the parameter range by devel- 
oping a worst-case and average case objective for optimizing the precision. Focusing on the single parameter 
case, we show that the optimization problems are linear programs. For the average case the solution to the linear 
program can be expressed analytically and involves a simple search: finding the largest element in a list. An 
example is presented which compares what is possible under constraints against the ideal with no constraints, 
the Quantum Fisher Information. 



INTRODUCTION 

The theoretical limit on the accuracy of parameter estima- 
tion in quantum metrology applications has been examined in 
depth, e.g., [1-6]. These studies reveal that special preparation 
of the instrumentation - the probe - can achieve an asymp- 
totic variance smaller than the Cramer-Rao lower bound [7], 
often referred to as the Quantum Fisher Information, abbrevi- 
ate here as QFI. In addition, the unique quantum property of 
entanglement can increase the parameter estimation conver- 
gence rate for N identical, independent experiments from the 
shot-noise limit of 1/ \/N to the Heisenberg limit of 1 /N. 

It is reasonable to expect, with or without entanglement, 
that the QFI will not be obtained with imperfect and limited 
instrumentation resources, i.e., not all states can be prepared 
and not all measurement schemes are possible. Under these 
conditions what exactly is the best that can be done? 

In this paper we present an approach which maximizes the 
parameter estimation accuracy in the presence of limits on in- 
strumentation, The method is based on the convex optimiza- 
tion approach to optimal experiment design as developed in 
[8] and as applied to quantum tomography in [9]. Incorpo- 
rating prior knowledge of the parameter range, we develop 
a worst-case and average case objective for optimizing the 
precision. Focusing on the single parameter case, we show 
that the optimization problems are linear programs. For the 
average case the solution to the linear program can be ex- 
pressed analytically and involves a simple search, i.e., find- 
ing the largest element in a list. This means that an enormous 
number of combinations of state and sensor configurations can 
be efficiently evaluated. 



OPTIMAL EXPERIMENT DESIGN 

Consider a quantum system dependent on an unknown 
scalar real parameter 9 which is known a priori to be in a 
set Q = { 9 | 9 min < 9 < 9 max }. The parameter 9 is to 
be estimated using data from repeated independent, identi- 
cal experiments. In each experiment the system can be put 
in any one of k = 1, . . . , iVconfig configurations. These 
are the available settings of input states and measurements. 
Each experiment in configuration k results in one of A^ut 



outcomes with probability P i\ k (9), i = l,---,Aout- Let 
Ni\k{9) denote the number of times outcome i is obtained 
from N k identical experiments in configuration k. Thus, 
Ei%(0) = N kPilk (6), EtT N i\k(0) = N k where E is 
the expected value operator with respect to the probability 
distribution Pi\k{9). Let A~ denote the total number of experi- 
ments and Afe the distribution of experiments in configuration 
k. Thus, A fe = N k /N => J2 k =T e A fc = 1 T A = 1. The 
problem is to select the distribution of experiments per config- 
uration, Afc, k = 1, . . . , A^config, or equivalently the number of 
experiments per configuration, N k , so as to obtain an estimate 
of 9 E with the best accuracy from N experiments. The 
"best" attainable estimation accuracy is defined here as the 
smallest possible Cramer-Rao bound on the estimation vari- 
ance [7]. 

Specifically, if 9n is an unbiased estimate of 9 from A" data, 
then the estimation error variance satisfies, 

NF(X, 6) E(9 N ~9) 2 > 1 

F(A,0) = A^) = E^T tig 2 ^, 9fc (0) (1) 

fflfe(0)=E£=T (VeP Ak {6)) 2 /p Ak {9) 

To achieve the best accuracy we will select A so as to maxi- 
mize a measure of the size of the Fisher Information, F(X, 9). 
To account for the knowledge that 9 £ we will consider two 
experiment design objectives for selecting A: average case 
and worst-case. 

Average-Case Experiment Design 

maximize F ac (X) = A T g avg (2) 
subject to 1 T A = 1, A^A is a vector of integers 

with g avg = J p(9)g(9)d9 where p(9) is the probability den- 
sity associated with 9 G 0. Although the objective function 
(average Fisher information) is linear in A, the integer con- 
straint on A makes the optimization problem hard. Utilizing 
the optimal experiment design method presented in [8, §7.5], 
the integer constraint is relaxed to the linear inequality A > 0. 
In addition, suppose we take a finite number of samples from 
the set 0, say, { 9 r \ r = 1, . . . , Ng } Then the non-convex in- 
teger optimization (2) is approximated by, 

maximize F ac (A) = X T g avg , g avg = ^ r p(&r)g{0 r ) m 
subject to 1 T A = 1, A > ( ' 
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FIG. 1: Quantum system for estimating parameter 9. 



This is a convex optimization problem in A, in fact, it is a 
linear program (LP). However, a particular advantage of this 
formulation (3), is that the solution is given explicitly by, 



A* 



1 k = argmaxfc' J2 r P(®r)gk>{9 r ) 
otherwise 



(4) 



with the optimal objective F ac (A) = maxfc ^ r p(9 r )gk{9 r ). 
It is possible that there is more than one optimal distribution 
because maxj may not be unique. However, due to limits 
on numerical precision, it is more likely that there are other 
choices which give similar results to the optimal objective. 

Worst-Case Experiment Design 

maximize F WC (A) = min^e ^ T g(9) (5) 
subject to 1 T A = 1, A^A is a vector of integers 

As in the average-case, relaxing the integer constraint and ap- 
proximating the objective function over a set of 9 sampled 
from the known set gives the optimization problem: 



maximize 



F wc(A) - 

subject to 1 T A = 1 



min r X T g( 
A > 



) 



(6) 



This is also an LP in A, but unlike the average-case, there is 
no explicit solution. However, it can be solved efficiently for 
a very large number of configurations N con fi g . A potential ad- 
vantage of the average-case solution over the worst-case solu- 
tion is that only a single configuration is required. As we will 
see in the example to follow, the two distributions can be quite 
different even though the Fisher information is similar. 

The solution to both of the relaxed and approximated prob- 
lems (3), (6) provide upper and lower bounds to the unknown 
solution of each with the integer constraint active. Specifi- 
cally, let A opt denote a solution to either (3) or (6) with the 
integer constraint. Let A be a solution to the relaxed (LP) 
versions. From the latter we can determine a nearby solu- 
tion which satisfies the integer constraint, e.g., set A = 
round(A). Then, F ac (\) < F ac (A opt ) < F ac (A rnd ) and 
Nk = NX" ld is the number of experiments to repeat in con- 
figuration k. 



Suppose that Q(9) can be described in terms of the Kraus Op- 
erator Sum Representation (OSR) with elements Qk{9). Then 
the outcome probabilities are: 

Pi{<j>,P,6) = Tr Mi{(j))(T{6, P), i=l,..., JV out 

<r%0) = Z k Qk(O)p(0)Qk{9y U} 

The state a(9, (3) is the output of the quantum chan- 
nel Q(6) and the input to the POVM. Suppose that 
the input and POVM configuration parameters can be 
selected, respectively, from { Pe \ I = 1, . . . , Mnput } and 
{ 4>u | k — 1, ... , N povm }. Hence, under the stated condi- 
tions, the worst-case experiment design problem (6) becomes, 



maximize min Y^k^i™ LfcT"' ffOfc, Pi, 9 r ) 

subject to X H > 0, Efrr EfcT *w = 1 

g(4>, P, 9) = EtT (Ve P l {4>, P, 0)) 2 /Pi{<j>, p, 9) 

(8) 

Similarly, the average-case experiment design problem (3) be- 
comes, 

maximize J2k=7 m EfcT"' ^ SWgO/'fe, Pi) 

subject to \ M > 0, Ef=°r EfcT x ki = 1 (9) 

SWgOfc,/^) = Y,rllP(®r)g(<i>k,Pe,9r) 

The worst-case distribution, A wc , is obtained by solving the 
LP (8). Following (4), the average-case distribution, A ac , 
which solves (9) is explicitly, 



\ ac 

A k£ — 



1 k,£ ~ argmaxg avg (0 fc /,/3 £ ') 

k',V (10) 

otherwise 



Solutions to (8) and (9), respectively, A wc and A ac , can be used 
to evaluate the worst-case and average-case levels of Fisher 
information as a function of the uncertain parameter 9 £ 0: 

f(x^,9) = EfrrE& ut >®9(<f>k,Pi,o) (id 

F(A ac ,fl) = Ajg g(<j)k,Pi,&) (12) 

In addition, as benchmarks for each 9 S 0, we can compute 
the maximum possible, subject to the constraints on the input 
and measurement scheme, and the QFI which is the maximum 
possible with no measurement constraints: the POVMs do not 
depend upon a configuration parameter as in (7). The maxi- 
mum subject to the constraints is, 



,(6) = max g{4> k ,Pt,9) 

h.t 



(13) 



QUANTUM SYSTEM PARAMETER ESTIMATION 

For the quantum system depicted in Figure 1, the quantum 
channel, Q(9), depends on the parameter 9 6 0, the input 
state, p(P), is dependent on the input configuration parame- 
ter P, and the POVM elements, Mj(0), i = 1, . . . , A^ut with 
E^Mi^) = /, depend on the configuration parameter <j>. 



For the single parameter system of Figure 1, the QFI is given 
by, [1,2], 



F QF1 (6,P)=TrS(6,p) 2 a(6,P) 

S(9, P)a(9, P) + (7(0, P)S{9, P) = 2V e a(9, P) 



(14) 



with a(9,P) from (7) and S(9,P) the solution to the above 
(matrix) Lyapunov equation. Fqyi(9, P), generally depends 
on the unknown parameter value 9, and in this case also on 
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the input configuration parameter (3. As developed in [1-6], 
for a unitary channel of the form U(9) = exp(—i9Ho), there 
is a ^-dependent pure state input \ij){ff)) such that the QFI is 
explicitly, 



F Q m(6) = (A 

max (H ) - Xmin(H )Y 



(15) 



with A max , A m j n here denoting the maximum and minimum 
eigenvalues of the Hamiltonian Hq. 

We ought to mention that the form of the system shown in 
Figure 1 is not the most general. For example, the "OSR" 
block might depend jointly on both 9 and a configuration pa- 
rameter a. The method, however, remains the same. 



EXAMPLE: PERTURBED UNITARY CHANNEL 

To illustrate the optimization methods we assume the quan- 
tum channel in Figure 1 is a unitary channel whose out- 
put is corrupted by amplitude damping. The unitary part 
is U{6) = exp(-i6H ), with H = -±= [ \ \ ] and with 
the unknown parameter 9 uniformly distributed in the set, 
6 = { 9 | 0.2 < (9/(7r/2) < 0.8 }. The amplitude damping 
channel can be described by an OSR with two elements (see, 
e.g., [10]), A x ( 7 ) = [ I A 2 (i) = [° ^] with 7 

the probability of dissipation. It follows that the OSR of Q(9) 
in Figure 1 has two elements, Qk{9) = Ak("f)U(9), k = 1,2. 

The available input for the experiment is the 2x1 pure state 
\ip((3)) which can be adjusted via an angle j3 as: \ip(f3)) = 
cos/3|0)+sin/3|l), < < tt. The POVMs can be adjusted 
via an angle <f> as: 



M 2 (<rj) 

k(0)> 



h-Mi{ft>) 

cos0|O) + sine 



< 6 < tt 



|1> 



We determine the Fisher information for two amplitude damp- 
ing probabilities: 7 £ {0, 0.25} with Ng = 100 uniformly 
spaced samples of 9 G 6. The POVM and input configura- 
tion angles f3, cf> are selected from their allowable ranges with 
-/Vinput = 10 and iVp 0vm = 10 uniformly spaced samples for 
each of the following three configuration constraints: 

1. POVM configured (0 < <j> < tt), input fixed (f3 = 0) 

2. POVM fixed (0 = 0), input configured (0 < (3 < tt) 

3. POVM & input configured (0 < /3 < tt, < 4> < tt) 

Figures 2-3 show the Fisher information as a function of the 
parameter 9 for the two values of amplitude damping and the 
three configuration constraints. In each figure the dotted lines 
are the QFI for each 9 (14). Note that the absolute maximum 
for the QFI is achieved only for 7 = (unitary channel) and 
using (15) with Hq as given above gives Fqfi(9) = 4. The 
solid lines are the maximum achievable for each value of 9 
that maximizes the Fisher information under the configuration 
constraints (13). The dashed lines are what is achieved by 



POVM configured, Input fixed, y = 



4r 
3 - 



1 — 
0.2 



0.3 0.4 0.5 0.6 0.7 

POVM fixed, Input configured, y = 



0.3 0.4 0.5 0.6 0.7 0.8 

POVM & Input configured, y = 



QFI at each e 

- max at each e 

- worst-case over 6 
average-case over 9 



0.5 
9/(jt/2) 



0.6 



FIG. 2: Comparison of configuration constraints with 7 = 
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FIG. 3: Comparison of configuration constraints with 7 = 0.25 



using the worst-case distribution of experiments (11), and the 
dot-dash lines are the average-case distribution of experiments 
(12). 

In all cases, the constrained Fisher information F(A ac , 9), 
and F(A WC , 9) are relatively close, sometimes nearly coinci- 
dent to the maximum possible, F max (9), and all are lower than 
the QFI. When both POVM and input are jointly configured 
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TABLE I: Optimal distributions 




FIG. 4: QFI vs. input configuration parameter j3 for 7 £ {0, 0.25}. 
A are the iVinput = 10 available. 

the constrained information begins to approach F max (9). The 
curves for the case where only the POVM is configured are 
generally below those where only the input is configured. 

The numerically non-zero elements of the worst-case and 
average-case optimal distributions for all the cases are shown 
in Table I. By construction, only one input configuration is re- 
quired for the average-case distribution (10). The worst-case 
distribution requires up to 3 configurations when 7 = 0.25. 
In this example the configuration angles remain relatively un- 
changed exhibiting some robustness to the amplitude damp- 
ing probability 7. The worst-case distributions change more 
significantly. Given the relatively close levels of Fisher in- 
formation for 9 £ O, it would seem more prudent to use the 
single-setting of input and POVM obtained from the average- 
case optimization. In most experiments there is a penalty in 
terms of time to reset the configurations. 

In this example the available input configurations effect the 
QFI. The solid lines in Figure 4 show the QFI for each value 
of 7 vs. the input configuration angle j3 for a large number 



of samples in the range. The QFI in this case is independent 
of 9. The triangles show the Mnput = 10 available values. 
The solid lines indicate that multiple inputs can achieve the 
bound whereas the restricted set forces a unique maximum 
which does not necessarily occur at the true maximum. For 
example, as seen in the top plot for 7 = 0, the constrained 
maximum is near the global maximum (Fqfj(9) = 4). This 
is achieved only in the case with 7 = and clearly over 
bounds the plot for 7 = 0.25. As might be expected, a per- 
turbation of the unitary channel, in this case via amplitude 
damping, makes it harder to attain the maximum possible QFI. 
Observe also that if the inputs were further constrained, say 
/3/tt e {0, 0.2, 0.5, 0.8}, then the achieved QFI would not be 
nearly as close to the maximum possible. The analysis of this 
examples thus provides the designer with information about 
the limit of performance of the system. If the potential perfor- 
mance increase over what is available under the constraints on 
instrumentation is significant, then a more flexible instrumen- 
tation might be considered worthwhile. 

CONCLUSION 

We have shown that maximizing the precision in estimating 
a single parameter in a quantum system subject to input and 
POVM constraints reduces to a linear program for both what 
is defined here as a worst-case and average-case objective. For 
the average-case, the solution to the linear program can be ex- 
pressed analytically and involves a simple search, i.e., find the 
largest element of an easily computed vector. Both solutions 
provide different levels of Fisher information over the range of 
anticipated parameter variation. Comparing these constrained 
solutions to the best possible under the constraints as well as 
to the QFI gives an indication of the performance limitations 
imposed by the constraints. 

Future efforts will consider the effect of entanglement and 
multi -parameter estimation. 
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